Quasiparticle interference in antiferromagnetic parent compounds of Fe-based 

superconductors. 
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Recently reported quasiparticle interference imaging in underdoped Ca(Fei_iCoa;)2 As2 shows 
pronounced C2 asymmetry that is interpreted as an indication of an electronic nematic phase with 
a unidirectional electron band, dispersive predominantly along the b-axis of this orthorhombic ma- 
terial. On the other hand, even more recent transport measurements on untwinned samples show 
near isotropy of the resistivity in the ab plane, with slightly larger conductivity along a (and not 6). 
We show that in fact both sets of data are consistent with the calculated ab initio Fermi surfaces, 
which has a decisevly broken C4, and yet similar Fermi velocity in both directions. This recon- 
ciles completely the apparent contradiction between the conclusions of the STM and the transport 
experiments. 

PACS numbers: 74.20. Pq,74.25.Jb,74.70.Xa 
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The Fe-based superconductors present a new paradigm 
for high-T(7 superconductivity as here Cooper-pairs ap- 
pear to emerge upon chemical doping from a metalhc 
ground state as opposed from a Mott insulator as found 
in the celebrated High-Tc cuprates^,. Despite this dif- 
ference of parent ground state of the Fe- and Cu-based 
superconductors, similarities he in that in both cases su- 
perconductivity emerges after the suppression of static 
ordered magnetism^. Although band theory has cor- 
rectly predicted the unusual antiferromagnetic (AFM) 
order in the parent compounds of the Fe-based super- 
conductors, it consistently overestimates the tendency to 
magnetism and underestimates the electronic mass, so 
there is no doubt that electronic interactions can not be 
ignored in quantitative descriptions, and that they play 
a different role compared to cuprates. The exact role of 
correlations, especially once the parent phase of the Fe- 
superconductors is doped, has been the focus of much 
debate and controversy. 

An almost universal feature of the Fe-supercondnctors 
is that in the parent phases, there is a tetragonal to 
orthorhombic structural phase transition that is closely 
associated with the onset of antiferromagnetic order—. 
Upon chemical doping x, the onset of the structural and 
magnetic transitions (T5 and T/v respectively) decrease 
with X and superconductivity emerges. The physical na- 
ture of the cross over from antiferromagnetic order to 
superconductivity varies between specific materials. In 
some cases both Ts and T/v coincide while in others T5 
is a few degrees higher than Tn^. 

Band structure calculations have suggested that the 
AFM ordering is accompanied by a strong restructuring 
of the Fermi surface, with the Fermi surface area being 
reduced by roughly an order of magnitude. This has been 
confirmed by optical and Hall measurements that regis- 
ter a drastic reduction of the carrier concentration in the 
AFM stated. The calculated AFM Fermi surface consists 
of several small pockets, which are arranged in the Bril- 
louin zone in a way that strongly breaks the tetragonal 
symmetry, but each of them is rather isotropio^. This 



led to a prediction of small transport anisotropy. An al- 
ternative point of view, that associates the orthorhombic 
transition with orbital (charge) degrees of freedom, sug- 
gests a double exchange (metallic) ferromagnetic interac- 
tion along one crystallographic direction and a superex- 
change along the other direction. This picture is also 
consistent with the observed AFM order and naturally 
suggests a metallic conductivity along the ferromagnetic 
chains and a substantially reduced conductivity in the 
other direction. 

Recent experiments on detwinned single crystals sup- 
port the former point of view: they demonstrate a small 
anisotropy with the AFM direction being more, not less 
metallic. However, transport measnrments are integrated 
probes, and also involve possibly anisotropic scattering 
rate, therefore experiments directly probing the topol- 
ogy of the Fermi surface in the AFM state are highly 
desirable. 

One such experiment has been recently performed by 
Chuang et al} . They have reported quasiparticle inter- 
ference (QPI) imaging of a lightly cobalt doped sample 
of CaFe2As2 compound. They interpreted their result in 
terms of a quasi-lD ("unidirectional") electronic struc- 
ture, metallic only along the FM, consistent with above- 
mentioned orbital picture. On the the other hand, their 
argumentation was rather indirect, based largly on the 
fact that directly measured dispersion of the QPI max- 
ima (which was indeed ID) coincded with the ARPES- 
measured band dispersion along the the same direction. 

In this paper we show that in reality the data of Ref. 
[l| are consistent with the calculated ab initio Fermi sur- 
faces, and not with the implied in that work ID bands. 
This reconciles completely the apparent contradiction be- 
tween the conclusions of Ref. 'r and the transport mea- 
surements on untwinned samples. 

The reported STM examination shows a QPI pattern 
in the momentum space that breaks completely the C4 
symmetry, the main features being two bright spots along 
the y (crystallographic 6) direction, with no counterparts 
along X (note that y is the ferromagnetic direction, and 



X in the antiferromagnetic one). Ref. [l| insists "that 
the scattering interference modulations are strongly uni- 
directional, which should occur if the k-space band sup- 
porting them is nematic" . However it should be kept in 
mind that this occurs in that part of the phase diagram 
where the long-range antiferromagnetic order is fully es- 
tablished, as reflected by the fact that the lattice sym- 
metry is orthorliombic, and the C2 symmetry is already 
completely broken. Indeed the size of the orthorhombic 
distortion is not "minute", as Ref. [l] posits, with b/a 
~ 1%, and is instead comparable with distortions seen in 
various iron oxides systems. For instance, in the Verwey 
transition the Fe-0 bond dilation is ~0.6% with Fe atoms 
in the same tetrahedral symmetry as in the ferropnictide 
superconductors'", and this is usually considered to be 
a strong distortion. Similarly, in the antiferromagnetic 
phase of FeO, where the cubic symmetry is completely 
broken, the structural effect is also on the same order™. 

Since the sample under study is orthorhombic it is mis- 
leading to call its electronic structure nematic, as the lat- 
tice orthorhombic distortion here is substantial. Nematic 
phases are frequently found in organic matter. The defin- 
ing characteristic of these phases is orientational order in 
the absence of long range positional order, resulting in 
distinctive uniaxial physical properties. It has also been 
proposed that nematic order exists in some electronic sys- 
tems, and may even play a role in mediating high tem- 
perature superconductivity^. Borzi et al^ demonstrated 
the presence of another interesting phase in Sr3Ru2 07 at 
millikelvin temperatures and high magnetic fields, which 
has also been called nematic. In this case, the crystallo- 
graphic planes were shown to remain strictly tetragonal 
(withing 0.01%) with C4 structural symmetry, while a 
pronounced C2 asymmetry in electronic properties was 
measured. This breaking of the electronic symmetry 
compared to that of the underlying lattice is now con- 
ventionally referred to as electronic nematicity (in fact, 
even in those cases one has to be careful to distinguish 
between nematic physics and simply an unusually weak 
electron-lattice coupling, but this goes beyond the scope 
of this paper, and in any event is not a concern for Fe 
pnictides where this coupling is strong). 

Since the tetragonal symmetry is decisively broken at 
the onset of the magnetic order in this ferropnictide, it 
is clear that the symmetry of the electronic structure 
defining the structural distortion is also completely bro- 
ken. What is more important is that while the observed 
QPI pattern does violate the C4 symmetry, it is clearly 
not one-dimensional, in the sense that it varies equally 
strongly along and ky directions. Thus, interpretation 
of the data in terms of a ID electron band does not ap- 
pear to be possible. To understand this experiment one 
needs to start with a realistic model for the electronic 
structure and actually calculate the QPI pattern. 

Such calculation has recently been presented by Knolle 
et al^. They used a weak-coupling theory that inter- 
prets tha antiferromagnetic state as resulting from a spin- 
Peierls transition, with a correspondingly small magnetic 
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FIG. 1: (Color online) Calculated LDA Fermi surface for 
CaFe2As2 in the antiferromagnetic state. 



moment. Knolle et al have been able to describe quali- 
tatively the experimental data obtained by Chuang et al 
in the sense that their calculated QPI pattern strongly 
breaks the C2 symmetry, while the band dispersion, on 
average, remains fairly isotropic in plane. Note that one 
should not be looking for a quantitative interpretation, 
since the STM experiment in question did not detect any 
Ca atoms on the surface, so the sample surface is likely 
charged with up to 0.5 hole per Fe, and thus any bulk 
calculation can only be applied to this experiment in a 
qualitative way. Besides, it was recently shown^^ that 
Fe pnictide systems feature surface states quite different 
from the bulk that should undoubtedly affect the STM 
spectra. 

However, this result, as mentioned, has been obtained 
in a weak coupling limit, corresponding to small magne- 
tization, while in this system the ordered magnetic mo- 
ments are on the order of 1 /is, and local moments even 
larger—"—. Not surprisingly, their Fermi surface is rather 
far from that measure recently on untwinned samples by 
Wang et a^, while the LDA Fermi surface reproduces 
it quite wellii. Indeed, this is a known problem in the 
weak coupling approach: while being physically justified 
for the paramagnetic parts of the phase diagram, the Fe 
magnetism in the ordered phases is driven by the strong 
local Hund rule coupling, and not by the Fermi surface 
nesting, as assumed in the weak copling models. 

Therefore we have calculated the QPI images for anti- 
ferromagnetic CaFe2As2 entirely from first principles'^ , 
using the Local Density Approximation (LDA) magnetic 
moment (somewhat larger that the experimental mo- 
ment at zero doping). We used the standard linear 
augnimented plane wave method as implemented in the 
WIEN2k codel^. The corresponding Fermi surface is 
shown in Fig. 1. We see that the magnetism has a dras- 
tic effect on the Fermiology, and the resulting Fermi sir- 
faces are completely breaking the C4 symmetry. Apart 
from small quasi-2D tubular pockets, originating from 
Dirac cones, there is one hole pocket around Z (0,0,7r/c 
or 27r/a,0,0) and two electron pockets between Z and 
0,7r/6,7r/c. It is immediately obvious that the QP scat- 
tering between these pockets must exhibit strong inter- 
ference for scattering along 6, but not a. 
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FIG. 2: (Color online) Quasiparticle interference pattern (in 
arbitrary units) for zero bias and qz~0, calculated using the 
same electronic structure as in Fig. [T]and Eq. 1. 



Indeed, we have calculated the QPI function Z, using 
the known expression (Ref. 19, Eq. S9) 

.^-^Y.5{E-E^)5{E'-E^+^), (1) 

k 

where we assumed a constant inpurity scattering rate 
and a constant tunneling matrix elements. This approx- 
imation is sufficient for a qualitative or semiquantitative 
comparison. As explained above, given that the surface 
in the experiment in question was charged compared to 
the bulk, a quantitative comparison is meaningless. 

A calculated pattern (there is some dependence on 
and on E, but we are interested in the qualitative features 
only) are shown in Fig. 2. One can see iimediately^that, 
very similarly to the patterns obtained in Ref. |l|, two 
sharp maxima appear at q = 0,±^,0, where ^ ~ ir/Ab. 
The origin of these QPI features is obvious from the 
Fermi surface (Fig. 1). Note that these LD A calculations 
have no adjustable parameters, and yet are in excellent 
qualitative agreement with the QPI images. 

It is also worth noting that while the calculated 
Fermi surfaces completely break the tetragonal sym- 
metry, which is fully reflected in the QPI images, the 
individual pockets are very three-dimensional, so that 
the calculated conductivity is comparable for all three 
directions^. While experimentally there is up to a 20% 
a/b charge transport anisotropy^ close to tetragonal to 
orthorhombic phase boundary in CaFe2As2, it is much 
less than what would be predicted for a quasi ID elec- 
tronic band, and of the opposite sign22. 

It may be worth at this point to explain at some length 
while a quantitative comparison between a Fourier trans- 
form of a tunneling current map, and theoretical calcu- 
lations, whether ours or any other, is impossible at this 
stage. Quasiparticle interference, as discussed in many 
papers, manifests itself in tunneling in a very indirect 



way. In a sense, it is a multistage process. First, a de- 
fect existing near the metal surface, is sdreened by the 
conducting electrons. This creates Friedel oscillations in 
the real space. This oscillations are formed by all elec- 
trons (mostly those near the Fermi surface, but not only) . 
In a multiband system, it includes electrons originated 
from different atomic orbitals, such as xy, xz, yz, z^ and 
x^ — y^ ■ As is well known in the theory of tunelling, the 
rate at which electrons tunnel through vacuum depends 
drastically on their orbital symmetry, especially on their 
parity (see, e.g., Ref)^). Indeed tunelling through a wide 
barrier mainly proceeds through electrons with zero mo- 
mentum projection onto the interface plane (such elec- 
trons have to travel the shortest lengths in the subbar- 
rier regime) . If such electrons belong to an odd 2D rep- 
resentation (for d-electrons, all but z^, if z is the normal 
direction), the tunneling rate is suppressed. This effect is 
well known in spintronics, where it can drastically change 
the current spin polarization. On the other hand, for a 
thin barrier the tunneling conductance depends on the 
number of the conductivity channels, which is given by 
the density of states (DOS) times normal velocity. In 
both cases, it is not just the density of quasiparticles, as 
assumed in Eq. 1 (and in Ref>ii), but the DOS weighed 
by a strongly k-dependent, unknown function. 

Nothing is known about the nature of the scatter- 
ing centers, producing the above mentioned Friedel os- 
cillations. In this particular experiment they may be 
magnetic or nonmagnetic defects, twin domain bound- 
aries, antiphase domain boundaries, remaining surface 
Ca ions, and more. Some of these scatterers are strongly 
anisotropic by nature, others are strongly dependent on 
the orbital character. We have dropped the scatter- 
ing matrix elements completely form our consideration. 
KnoUe et al^ instead have chosen a specific model for the 
scattering centers. We believe that without any knowl- 
edge about the actual scattering centers in the system 
any QPI using a particular model is more obscuring the 
actual physics, compared to the simplest constant matrix 
elements approximation, rather than clarifying it. 

Finally, there are several issues specific for this par- 
ticular experiment: (1) unknown, but strongly different 
from the bulk, charge state. As opposed to Bal22, and 
Srl22, where 1/2 of the alkaline earth atoms stay on the 
surface, providing charge neutrality, in Cal22 STM does 
not detect any Ca on the surface, suggesting a strongly 
charged surface. A corollary of that is appearence of a 
surface reconstruction (as indeed observed), of a surface 
relaxation, and, importantly (since tunneling proceeds 
largely through the surface states) , of surface bands (as 
demonstrated, for instance, in Ref.^^. 

While the above considerations preclude a quantitative 
comparison and extracting quantitative analysis of the 
experiment in question, we see, particularly when com- 
paring our calculations with those of Knolle et ai^, that 
the C2 QPI structure observed in Ref. [l| is a very univer- 
sal consequence of the long-range stripe-type antiferro- 
magnetic ordering. Indeed, Knolle et al calculations were 
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built upon a besically incorrect band structure and fermi 
surfaces, an used a weak coupling nesting scenario for 
the antiferromagnetism, while in reality the magnetism in 
pnictides is a strong coupling phenomenos; yet, their cal- 
culations produced a "unidirectional" QPI pattern just 
as well. Together with the strong-coupling LDA calcu- 
lations, this span a large range of possible models, indi- 
cating that the C4 symmetry is strongly broken in QPI 



images with simply by virtue of the long range AFM or- 
der, whatever the the origin of this order. 

Last but not least, we can also predict, from our cal- 
culations, that this symmetry will be also broken, al- 
though the peaks are likely to be substantially broaden, 
in the truly nematic phase (see reviewi^ for a discussion), 
that is to say, the phase between the long-range magnetic 
transition and the structural orthorhombic transition. 
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